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Comparison of calculated and experimental isotope edited 
FTIR difference spectra for purple bacterial photosynthetic 
reaction centers with different quinones incorporated into 
the Qa binding site 



Previously we have shown that ONIOM type (QM/MM) calculations can be used to 
simulate isotope edited FTIR difference spectra for neutral ubiquinone in the Qa binding 
site in Rhodobacter sphaeroides photosynthetic reaction centers. Here we considerably 
extend upon this previous work by calculating isotope edited FTIR difference spectra 
for reaction centers with a variety of unlabeled and 18 0 labeled foreign quinones 
incorporated into the Qa binding site. Isotope edited spectra were calculated for reaction 
centers with 2,3-dimethoxy-5,6-dimethyl-1 ,4-benzoquinone (MQo), 2,3,5,6-tetramethyM , 
4-benzoquinone (duroquinone, DQ), and 2,3-dimethyl-l,4-naphthoquinone (DMNQ) 
incorporated, and compared to corresponding experimental spectra. The calculated and 
experimental spectra agree well, further demonstrating the utility and applicability of our 
ONIOM approach for calculating the vibrational properties of pigments in protein binding 
sites. The normal modes that contribute to the bands in the calculated spectra, their 
composition, frequency, and intensity, and how these quantities are modified upon 18 0 
labeling, are presented. This computed information leads to a new and more detailed 
understanding/interpretation of the experimental FTIR difference spectra. Hydrogen 
bonding to the carbonyl groups of the incorporated quinones is shown to be relatively 
weak. It is also shown that there is some asymmetry in hydrogen bonding, accounting 
for 10-13 cm -1 separation in the frequencies of the carbonyl vibrational modes of the 
incorporated quinones. The extent of asymmetry in H-bonding could only be established 
by considering the spectra for various types of quinones incorporated into the Qa binding 
site. The quinones listed above are "tail-less." Spectra were also calculated for reaction 
centers with corresponding "tail" containing quinones incorporated, and it is found that 
replacement of the quinone methyl group by a phytyl or prenyl chain does not alter ONIOM 
calculated spectra. 
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INTRODUCTION 

Quinones play an important role in biological proton and elec- 
tron transfer processes that occur in both respiration and photo- 
synthesis (Trumpower, 1982). In type II photosynthetic reaction 
centers two quinone molecules act as terminal electron accep- 
tors (Ke, 2001a,b). The two quinones are often termed Qa and 
Qb- In this manuscript we will refer to the quinone binding site 
as Qa and Qb, however. The quinones that occupy the Qa and 
Qb binding sites have very different functions. The Qa quinone 

Abbreviations: Ala, alanine; DS, difference spectra/spectrum spectroscopy or 
spectroscopic; C— O, carbonyl; DQ, duroquinone or 2,3,5,6-tetramethyl-l,4- 
benzoquinone; DFT, density functional theory; DMNQ, dimethyl naphtho- 
quinone; ET, electron transfer; FTIR, Fourier Transform Infra-red; H-bond, hydro- 
gen bond; His, histidine; MQo, 2, 3-dimethoxy-5,6-dimethyl-l, 4-benzoquinone; 
MM, molecular mechanics; NQ, naphthoquinone; PQ, plastoquinone; PBRC, pur- 
ple bacterial reaction center; PED, potential energy distribution; Rh, Rhodobacter; 
QM, quantum mechanics; UQ„, ubiquinone-n; VK, vitamin Ki (phylloquinone). 



is an intermediary cofactor involved in transferring electrons 
from (bacterio) pheophytin to Qb, while the Qb quinone couples 
proton and electron transfer processes (Ke, 2001b, c,d). 

In this manuscript we focus on the Qa binding site. The 
quinone occupying the Qa binding site is species dependent. 
In Rhodobacter (Rb.) sphaeroides purple bacterial reaction cen- 
ters (PBRCs) a ubiquinone (UQ) molecule occupies the Qa 
binding site. In PBRCs from Blastochloris viridis (Shopes and 
Wraight, 1985) and Chloroflexus aurantiacus (Hale et al., 1983) 
a menaquinone occupies the Qa binding site. In photosystem 
II reaction centers from oxygen evolving organisms, a plasto- 
quinone (PQ) molecule occupies the Qa binding site. In photo- 
system I the secondary electron acceptor, termed Ai, is a vitamin 
Ki (VK) molecule (also called phylloquinone). Menaquinone and 
VK are both naphthoquinone (NQ) moieties that differ only in 
the degree of saturation of the tail at C6- 
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Figure 1 shows the structure, numbering and abbreviations we 
will use for the various quinones discussed in this manuscript. 
MQo and DMNQ are UQ and VK analogues, respectively, in 
which the hydrocarbon tail has been replaced with a methyl 
group. DQ is a PQ analogue in which the hydrocarbon chain at 
C6 is replaced with a methyl group. 

It has been suggested, at least for the UQ that occupies the 
Qa site in PBRCs, that the role of the hydrocarbon chain at 
is to anchor and orient the quinone head-group in a specific way 
(Warncke et al., 1994). Data is available that may argue against 
this proposal, however, (McComb et al., 1990; Srinivasan and 
Golbeck, 2009; Wraight and Gunner, 2009). 

Comparison of the properties of PBRCs with MQo and UQ„, 
or VK and DMNQ, incorporated into the Qa binding site will 
allow one to assess how or if the hydrocarbon chain at modi- 
fies the quinones functional properties. Similarly, comparison of 
the properties of PBRCs with MQo and DQ incorporated into the 
Qa binding site will allow one to assess how or if the methoxy 
groups at C2 and C3 modifies the quinones functional proper- 
ties. Calculated spectra for Rb. sphaeroides PBRCs with VK in 
the Qa binding site can be compared to experimental spectra. 
These calculated spectra may also serve as a useful model for B. 
viridis PBRCs that naturally have VK incorporated into the Qa 
binding site. 

Fourier transform infrared (FTIR) difference spectroscopy 
(DS) is a sensitive molecular-level probe of pigment-protein inter- 
actions, and it is widely used to study both the neutral and 
reduced states of quinones in PBRCs (Breton and Nabedryk, 
1996) and in photosystem II (Noguchi and Berthomieu, 2005). 
In this manuscript we focus on Q^/Qa FTIR DS. Many molecu- 
lar species contribute to Q^/Qa FTIR DS, and in the past it has 
been difficult to identify which bands are associated specifically 
with UQ in the Qa site. However, fully functional quinones can be 
incorporated into Qa depleted PBRCs, and by collecting Q^/Qa 



FTIR DS using PBRCs with unlabeled and isotopically labeled 
quinones incorporated, so called isotope edited FTIR difference 
spectra can be constructed, and from these spectra it has proven 
possible to separate contributions of the quinones from those of 
the protein in Q A /Q A FTIR DS (Breton and Nabedryk, 1996). 
Previously, a variety of unlabeled and ls O labeled quinones have 
been incorporated into the Qa binding site in PBRCs, and ls O 
isotope edited FTIR DS have been obtained (Breton et al., 1994a). 
The goal in this manuscript is the simulation of these experi- 
mental ls O isotope edited FTIR DS associated with the neutral 
state of the quinone in the Qa binding site. Calculated IR spec- 
tra associated with the quinone anion state are considerably more 
complicated (Lamichhane and Hastings, 2013) and are currently 
being undertaken. 

Although experimental Q^/Qa FTIR DS have been obtained 
using PBRCs with various unlabeled and isotope labeled 
quinones incorporated, virtually no computational work aimed 
at modeling the experimental FTIR DS have been undertaken. 
Calculations aimed at modeling the vibrational properties of 
quinones in the Qa binding site must account for the protein 
environment. Quantum mechanical (QM) calculations of such 
a large molecular system (pigment plus protein environment) 
are unfeasible. To circumnavigate this problem methods have 
been developed that allow one to separate the molecular system 
into distinct layers that can be treated at different levels of the- 
ory. In one layer the chemical properties of the principle species 
of interest (the pigment) can be calculated using "high-level" 
QM methods. The surrounding protein environment is included 
in the calculation but it is treated using computationally less 
expensive molecular mechanics (MM) methods. 

Recently, we have undertaken QM:MM calculations for UQ 
in the Qa binding site using the ONIOM method (Vreven et al., 
2006). ONIOM is an acronym for: our Own N-layered Integrated 
molecular Orbital + Molecular mechanics package. In these cal- 
culations we showed that we could simulate experimental isotope 
edited FTIR difference spectra obtained using PBRCs with neutral 
UQ in the Qa binding site. Here we extend upon these previ- 
ous studies by attempting to simulate experimental isotope edited 
FTIR difference spectra obtained using PBRCs with symmetric 
tail-less and corresponding tail containing quinones incorporated 
into the Qa binding site. We show that the calculated spectra agree 
well with the experimental spectra, further supporting the notion 
that the ONIOM method is a useful approach for understand- 
ing complex FTIR difference spectra associated with pigments in 
protein binding sites. We are also able to assess to what extent the 
quinone hydrocarbon chain may influence the calculated spectra. 

MATERIALS AND METHODS 
STRUCTURAL MODEL USED IN CALCULATIONS 

The molecular model used in ONIOM calculations was gener- 
ated from the crystal structure of Rb. sphaeroides PBRCs at 2.2 A 
resolution (Stowell et al, 1997) (PDB file 1AIJ). From the PDB 
file all atoms within 10 A of either carbonyl oxygen atom of UQ 
were selected. This subset of atoms formed the basis of the Qa 
binding site structural model. Hydrogen atoms (not included 
in the PDB file) were added to the model using the software 
Gaussview4, resulting in a final structural model consisting of 



UQ n MQo DQ 




DMNQ VK t 



FIGURE 1 | Structure and numbering of ubiquinone (2,3-dimethoxy, 
5-methyl,6-prenyl benzoquinone) (UQ„), 2,3-dimethoxy, 5,6-methyl 
benzoquinone (MQo), 2, 3, 5,6-methyl benzoquinone (duroquinone, 
DQ), 2,3-dimethyl, 1,4-naphthoquinone (DMNQ) and 2-methyl, 
3-phytyl 1,4-naphthoquinone (VK). The numbering scheme employed 
here for the naphthoquinone structures is nonstandard, and was chosen to 
facilitate comparison between naphthoquinone and ubiquinone structures. 
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1024 atoms. Following the addition of hydrogen atoms the struc- 
tural model was optimized (energy minimized) using ONIOM 
methods with all atoms associated with the protein backbone, 
and the non-heme iron atom, being held fixed. All atoms of the 
amino acid side chains and the incorporated quinone are uncon- 
strained. For calculations of UQ/VK the hydrocarbon tail was 
modeled as a single prenyl/phytyl unit, respectively. Inclusion of 
further prenyl/phytyl units did not alter the calculated spectra 
(not shown). 

The structural models for the different incorporated quinones 
are initially set up by simply replacing the molecular substituents 
of the originally incorporated UQ species. So, the C=0 groups 
of the different quinones incorporated will initially have the same 
orientation and position to that found for UQ in the Qa binding 
site. DMNQ and VK structures were constructed starting from 
the UQ structure, by replacing the methoxy groups with the NQ 
aromatic ring, without alteration of the quinone ring (Figure 2F). 
This orientation of the NQ ring (of DMNQ and VK) was cho- 
sen as previous docking calculations have suggested it is the most 
favorable (Hucke et al., 2002). 

CALCULATIONS 

All calculations were undertaken using Gaussian 03 (Frisch et al., 
2004) software. For calculation of UQ molecules in the gas 
phase, and for the QM part of ONIOM calculations, molecular 
geometry optimizations and harmonic vibrational frequency cal- 
culations were undertaken using hybrid DFT methods, employ- 
ing the B3LYP functional and the 6-31+G(d) basis set. This 
choice of functional and basis set are appropriate for calcula- 
tion of the vibrational properties of quinones (Wheeler, 2001; 
Bandaranayake et al, 2006). The MM part of the ONIOM 
calculation is undertaken using the AMBER force field (Case 
et al., 2005). Following ONIOM geometry optimization of 
the structural model, the optimized quinone molecule from 
the model is considered separately for vibrational frequency 
calculations. 

NORMAL MODE ASSESSMENT 

Assignment of calculated vibrational frequencies to molecular 
groups is based on a consideration of the calculated atomic 
displacements (in Cartesian coordinates) associated with the nor- 
mal modes. These atomic displacements can be animated using 
software (Gauss View4), and the molecular groups that most 
prominently contribute to the normal modes can be assessed 
visually [see Figure 2 in Bandaranayake et al. (2006)]. In addi- 
tion, potential energy distributions (PEDs) of the normal modes 
are calculated using the freeware GAR2PED (Martin and Van 
Alsenoy, 1995). For the various quinones in the Qa binding site 
we calculate both vibrational mode frequencies and intensities. 
With both the frequency and intensity information "stick spec- 
tra" can be constructed. These stick spectra are representative of 
IR absorption spectra, as described previously (Bandaranayake 
et al., 2006; Parameswaran et al., 2008; Lamichhane et al., 2010). 
The calculated stick spectra are convolved with a Gaussian func- 
tion with half- width of 4 cm -1 to produce more realistic looking 
IR absorption spectra (Bandaranayake et al., 2006; Parameswaran 
et al., 2008; Lamichhane et al, 2010). 



RESULTS 

UQ STRUCTURE AND NUMBERING 

Figure 1 shows the structure and numbering scheme used here 
for UQ„, MQ 0 , DQ, DMNQ, and VK. NQ's generally have a 
different numbering scheme. We have applied the UQ number- 
ing scheme to DMNQ and VK for the sake of easy comparison. 
Figure 2 shows a picture of (A) DMNQ, (B) VK, (C) MQ 0 , and 
(D) UQi in the Qa binding site along with the two H-bonding 
amino acids. The structures shown are after geometry optimiza- 
tion using ONIOM methods. Possible hydrogen bonds (or ligand 
to the non-heme iron atom) are indicated by dotted lines. 

To gain a better sense of the relative orientation of the dif- 
ferent quinones in the Qa binding site Figure 2 also shows 
the (E) DMNQ/DQ and (F) DMNQ/UQ structures from two 
ONIOM calculations overlapped. These overlapped structures are 
created by considering the (fixed) atoms of the protein back- 
bone. Figures 2E,F indicates that the side chains of HisM219 and 
AlaM260 are unaltered when a different quinone is incorporated 
into the binding site. 

Table 1 lists several bond lengths and bond angles derived 
from our ONIOM calculated optimized geometries of the var- 
ious quinones in the Qa binding site. For comparison, Table 1 
also lists corresponding bond lengths and angles derived from 
our DFT calculated optimized geometries of the various quinones 
in the gas phase. Table 1 also list results obtained from previ- 
ous QM/MM calculations (Sinnecker et al., 2006), and data taken 
from the 1AIJ crystal structure (Stowell et al., 1997). 

Table 1 list the three distances associated with the peptide or 
imidazole H-bond to the carbonyl oxygen atoms of the quinone 
(N-H, NH-O and N-O distances). These distances completely 
determine the N-H-O H-bond geometry. Associated angles are 
also listed in Table 1. 

For UQio in the Qa site the 1AIJ crystal structure (Stowell 
et al., 1997) indicates that the Ci=0 bond is marginally longer 
than the Gi=0 bond. (1.234 vs. 1.232 A). In contrast, from 
ONIOM calculations of all the quinones listed in Table 1 the 
Ci=0 bond is shorter than the C4=0 bond. 

In gas phase calculations the Ci=0 and C4=0 bond lengths 
are shorter than that found in ONIOM calculations (except the 
Ci=0 of DQ). This lengthening of the C=0 bonds of the 
quinones in the Qa binding site is related to hydrogen bond- 
ing and other electrostatic interactions of the pigment with the 
protein environment. 

For UQi and MQ 0 the Ci=0 bond is shorter than the C 4 =0 
bond in both ONIOM and gas phase calculations. In gas phase 
calculations this difference in Ci=0 and C 4 =0 bond lengths of 
UQi and MQo must relate to the differing orientations of the 
C2 and C3 methoxy groups. The ONIOM calculated C2 and C3 
methoxy group dihedral angles for UQi are —25.3 and 150.5° 
(Table 1). Similar angles are calculated for MQo. The calculated 
dihedral angles for UQi in the gas phase, and the observed angles 
for UQ10 in the Qa binding site (from the crystal structure) are 
within 32° of that calculated for UQi using ONIOM methods. 

For UQ10 in the Qa site the crystal structure indicates that the 
C 2 =C 3 bondis shorter than the C 5 =C 6 bond. (1.404 vs. 1.419A). 
In contrast, for all quinones except DQ, in both ONIOM and gas 
phase calculations the C2=C3 bond is longer than the Cs=C6 
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FIGURE 2 | Ball and stick representation of the calculated structure of 
(A) DMNQ, (B) VK, (C) MQ 0 and (D) UQ t in the Q A binding site in 
PBRCs from 17b. sphaeroides. Structure shown is after geometry 
optimization using ONIOM methods. Possible H-bonds are shown (dotted). 
Hydrogen atoms, except the ones involved in H-bonding, have been 




omitted. Distances quoted are in A. (E,F) Overlap of structures obtained 
from ONIOM calculations for (E) DMNQ and DQ and (F) DMNQ and UQi. 
Overlap is based on using the (fixed) backbone atoms of HisM219 and 
AlaM260. The amino acid side chains are virtually unaltered in the 
structures shown. 



bond. For VK and DMNQ the calculated C 2 =C 3 and C 5 =C 6 
bond lengths are considerably different. 

For UQi and VK the hydrocarbon chain attached at C6 makes 
a distinct "kink" after the first carbon atom (Figures 2B,D). The 
C-C-C bond angle is 112-115° for both quinones in both the 
ONIOM and gas phase calculations. The calculated angles are vir- 
tually the same as that found in the crystal structure. Given these 
similarities (between the ONIOM and gas phase calculations for 
both UQi and VK, as well as between calculation and experiment) 
the suggestion is that the protein environment does not constrain 
the orientation of the quinone ring relative to the C6 hydrocarbon 
chain. 

Figure 3 A shows ONIOM calculated ls O isotope edited differ- 
ence spectra for neutral VK and DMNQ. Figure 3D shows corre- 
sponding DFT calculated spectra for VK and DMNQ in the gas 
phase. Experimental spectra are also shown (Figures 3B,C) for 
comparison. Positive/negative bands in the isotope edited spectra 
are due to the unlabeled/ 18 O labeled quinone species, respectively. 
The ONIOM calculated spectra clearly better describe the experi- 
mental spectra. Calculations including the protein environment 



appear to be necessary in order to adequately simulate the 
experimental spectra. 

In the DFT calculated ls O isotope edited spectrum for 
DMNQ/VK in the gas phase (Figure 3D) the antisymmet- 
ric vibration of both C=0 groups gives rise to the band 
at ~ 1654 cm -1 (Figure 3D), which downshifts 30 cm -1 upon 
ls O labeling (Bandaranayake et al., 2006). The calculated isotope- 
edited gas phase spectrum is in excellent agreement with the 
experimental spectrum for DMNQ in solution (Breton et al., 
1994a). 

Figure 3 shows that except for a small frequency shift in some 
of the modes, the calculated spectra for DMNQ and VK are vir- 
tually the same. Replacement of the methyl group at C6 with 
an isoprene unit therefore, has no influence on the calculated 
spectra. 

The normal modes (frequencies and intensities) that give rise 
to the various bands in the ONIOM calculated ls O isotope edited 
spectra of DMNQ and VK are listed in Table 2. The PEDs, which 
quantify to what extent various internal coordinates contribute to 
the normal modes, are also listed in Table 2. Results for DMNQ 



Frontiers in Plant Science [ Plant Physiology 



August 2013 | Volume 4 | Article 328 | 4 



Zhao et al. 



Calculated Q a /Qa FTIR difference spectra 



Table 1 | Comparison of bond lengths and angles derived from the ONIOM calculated (0) and gas phase calculated (GP) optimized geometry of 
neutral UQ^ MQ 0 , DQ, VK, and DMNQ. 
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Bone/ lengths and angles from the 1AIJ crystal structure (Stowell et al., 1997) are also listed. Distances are in A and angles are in degrees. The C2 and C3 methoxy 
group dihedral angles are defined as the C3-C2-O-CH3 and C2-C3-O-CH3 dihedral angles from Sinnecker et al. (2006). 
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FIGURE 3 I (A) ONIOM calculated 1s O isotope edited DDS for neutral VK 
(dotted) and DMNQ (solid) in the Qa binding site. Experimental spectra 
for VK (B) and DMNQ (C) are also shown, and were taken from Breton 
et al. (1994a), with permission. (D) DFT calculated 18 0 isotope edited 
DDS for neutral VK (dotted) and DMNQ (solid) are also shown. 
ONIOM/gas phase calculated spectra were scaled by 0.9718/0.9608, 
respectively. 



and VK are very similar. Below we will discuss calculated data 
obtained for DMNQ with the recognition that very similar results 
and conclusions also apply to VK. 

In the ONIOM calculated ls O isotope edited spectra for 
DMNQ the two bands at 1654 and 1644 cm -1 are due to Q=0 
and C4=0 stretching vibrations, respectively (Table 2). These 
bands almost certainly correspond to the at 1653 and 1643 cm -1 
bands in the experimental spectrum (Figure 3C). So our cal- 
culations predict that the 1653 and 1643 cm -1 bands in the 
experimental spectrum are due to Q=0 and C4=0 stretching 
vibrations, respectively. This is in fact the first direct evidence that 
the 1653 and 1643 cm -1 bands in the experimental spectrum are 
due to Ci=0 and C4=0 stretching vibrations. A direct assign- 
ment has never been made because specific 13 Q and 13 C4 isotopic 
labeled VK or DMNQ has never been incorporated into the Qa 
binding site in PBRCs. 

In contrast to the observation of two separate C=0 
modes in ONIOM calculations, in gas phase calculations the 
two C=0 modes of DMNQ are anti-symmetrically coupled 
(Bandaranayake et al, 2006), giving rise to a single intense band 
at 1654 cm -1 (Figure 3D). 

In ONIOM calculations it is found that upon ls O labeling the 
Ci=0 and C 4 =0 modes of DMNQ downshift 23 and 19 cm" 1 , 
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Table 2 | Normal mode frequencies (in cm 1 ), intensities (in km/mol) and PEDs (in %) calculated using ONIOM methods for unlabeled and 18 0 
labeled neutral DMNQ, VK, DQ, MQ 0 , and UQ^ 
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Frequency shifts upon 18 O labeling are also listed. Negative signs in the PEDs refer to the relative phase of vibration of the internal coordinates. Only internal 
coordinates that contribute at least 5 % are shown. Mode frequencies were scaled by 0.9718. 



to 1631 and 1625 cm -1 , respectively (Table2). Upon ls O label- 
ing the 1625 cm -1 mode is more than six times more intense 
than the 1631 cm -1 mode. The 1625 cm -1 mode in ls O labeled 
DMNQ is due mainly to the antisymmetric coupled vibration of 
both C=0 groups. That is, two separate C=0 modes of unlabeled 
DMNQ couple upon ls O labeling. This behavior is not predicted 
based upon consideration of the experimental spectra, where it is 
"assumed" that the two C=0 modes remain separate upon ls O 
labeling (Breton et al., 1994b; Breton, 1997). 

A C=C mode of the quinone ring of DMNQ is found at 
1617 cm" 1 . This quinonic C=C mode downshifts 13 cm , to 
1604 cm -1 , upon ls O labeling, with little change in intensity. 
The 1604 cm -1 mode composition in ls O labeled DMNQ dis- 
plays considerable mixing of the C=0 and C=C modes (the 
C=0 modes account for 39% of the PED). This behavior might 
be expected given that the C=0 and C=C modes are closer in 
frequency upon ls O labeling. 

A C=C mode of the aromatic ring of DMNQ occurs at 
1591 cm -1 . This aromatic C=C mode downshifts 5cm _1 to 
1586 cm -1 upon ls O labeling. The 1586 cm -1 mode in ls O 



labeled DMNQ displays some mixing with C=0 modes (23%), 
and the intensity of the aromatic C=C mode nearly doubles upon 
ls O labeling. The calculated C=C normal modes and their inter- 
pretation in terms of internal coordinates, as well as the calculated 
18 O induced frequency shifts are similar to that suggested on the 
basis of the experimental spectra (Breton et al., 1994a). 

Figure 4A shows ONIOM calculated 18 0 isotope edited IR dif- 
ference spectra for neutral DQ in the Qa binding site. Figure 4C 
shows the corresponding DFT calculated spectrum for DQ in the 
gas phase. The experimental spectrum is shown in Figure 4B. 
Again, the ONIOM calculated spectrum agrees well with the 
experimental spectrum while the calculated gas phase spectrum 
does not. The ONIOM calculated normal modes (frequencies and 
intensities) that give rise to the bands in the ls O isotope edited 
spectrum for DQ, as well as the PEDs, are listed in Table 2. 

In the ONIOM calculated ls O isotope edited spectrum for 
DQ the two bands at 1646 and 1632 cm -1 are due to Ci=0 
and C4=0 stretching vibrations, respectively (Table 2). Upon ls O 
labeling the C 4 =0 and Q=0 modes both downshift 33 cm -1 
with little change in mode intensities (Table 2). This 33 cm -1 
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FIGURE 4 | (A) ONIOM calculated 18 0 isotope edited DDS for neutral DQ 
in the Qa binding site. Experimental spectra are shown in (B), and were 
taken from Breton et al. (1994a) with permission. (C) DFT calculated 18 0 
isotope edited DDS for neutral DQ. ONIOM and gas phase calculated 
spectra were scaled by 0.9718 and 0.9608, respectively. 



downshift is large compared to that calculated for DMNQ (19- 
23 cm -1 ). For DQ, the ONIOM calculated mode composition is 
virtually unchanged upon ls O labeling. This is also markedly dif- 
ferent to that calculated for DMNQ. The normal modes that give 
rise to the isotope edited spectra of DQ and DMNQ in the gas 
phase are very similar, so replacing the methyl groups of DQ with 
the aromatic ring of DMNQ in the Qa binding site leads to an 
alteration in the electronic structure of the quinone ring. This 
modification is not obvious given the similar orientation of DQ 
and DMNQ in the Qa binding site (Figure 2E). 

In the DFT calculated ls O isotope edited spectrum for DQ in 
the gas phase, the two C=0 modes are strongly coupled, and give 
rise to the band at 1639 cm -1 in Figure 4C. The calculated gas 
phase spectrum is in line with the experimental FTIR spectrum 
for DQ in solution (Breton et al, 1994a). 

From ONIOM calculations for DQ, C=C modes do not cou- 
ple with the C=0 modes. The C2=C3 and Cs=C6 groups couple 
to produce in phase and out of phase vibrational modes. The in 
phase mode has negligible IR intensity. The out of phase C=C 
mode is calculated to be at 1620 cm -1 . This C=C mode is vir- 
tually unaltered in frequency, intensity, and mode composition 
upon ls O labeling (Table 2). 

Figure 5A shows ONIOM calculated ls O isotope edited IR 
spectra for neutral MQn (solid) and UQi (dotted) in the Qa 
binding site. Figure 5D shows DFT calculated spectra for MQo 
and UQi in the gas phase. The experimental spectra for UQi 
and MQo are shown in Figures 5B,C, respectively. The normal 
modes (frequencies and intensities) that give rise to the various 
bands in the ONIOM calculated spectra, as well as the PEDs, 
are listed in Table 2. The data for UQi has been presented pre- 
viously (Lamichhane and Hastings, 2011), and we show here 
that very similar spectra are calculated for both MQo and UQi. 
Replacement of an isoprene unit at with a methyl group does 
not greatly alter the calculated spectra. 

In the ONIOM calculated ls O isotope edited spectra for neu- 
tral MQo the two bands at 1666 and 1627 cm -1 (Figure 5C) 
are due to Q=0 and C4=0 stretching vibrations, respectively 




Wavenumber (cm 1 ) 



FIGURE 5 | (A) ONIOM calculated ls O isotope edited DDS for neutral MQ 0 
(solid) and UQ q (dotted) in the Q A binding site. Experimental spectrum for 
(B) MQo and (C) UQ! are also shown, and were taken from Breton et al. 
(1994a) with permission. (D) DFT calculated 18 0 isotope edited DDS for 
neutral MQ 0 (solid) and UQ-| (dotted). ONIOM and gas phase calculated 
spectra were scaled by 0.9718 and 0.9608, respectively. 



(Table 2). Upon ls O labeling the Q=0 and C4=0 modes down- 
shift 35 and 40 cm -1 , respectively (Table 2). The mode intensities 
and composition are little altered by ls O labeling. The C=C 
mode of MQo/UQi at 1601 cm -1 up-shifts 8/7 cm -1 upon ls O 
labeling (Table 2). An explanation for this ls O induced fre- 
quency upshift has been presented (Lamichhane and Hastings, 
2011). 

In the DFT calculated ls O isotope edited spectra for UQi/MQo 
in the gas phase, the two C=0 modes are well separated, and 
give rise to bands at ~1665 and ~1635 cm -1 (Figure 5D), which 
downshift 36 and 27 cm -1 , respectively, upon ls O labeling, as 
described previously (Lamichhane et al., 2011). 

DISCUSSION 

Previously we have shown that ONIOM methods can be used 
to calculate isotope edited difference spectra for UQi in the Qa 
binding site, and that these calculated spectra model very well the 
corresponding experimental spectra (Lamichhane and Hastings, 
2011). Here we considerably extend these studies, and show that 
ONIOM calculated isotope edited spectra for different quinones 
in the Qa binding site model very well the corresponding exper- 
imental spectra. We further show that the calculated spectra for 
quinones in the gas phase are totally inappropriate for modeling 
the vibrational properties of quinones in the Qa binding site. This 
is also very likely to be the case for modeling the properties of any 
protein bound pigment. 

Without normal mode vibrational frequency calculations 
interpretation of experimental spectra is limited, and here we 
clearly show that our computational methods lead to a greatly 
increased understanding of the normal modes that contribute to 
the bands in the experimental spectra. 

In our ONIOM calculations for neutral UQi in the Qa bind- 
ing site we considered all amino acids included in the model 
only at the molecular mechanics level of computation. Previously, 
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ONIOM calculations have been undertaken in order to model 
EPR data associated with the UQ anion in the Qa binding site 
(Lin and O'Malley, 2008). In these calculations key amino acids, 
such as HisM219 and AlaM260, were considered at the higher 
quantum mechanical level of calculation. Taken together these 
results may indicate differences in the nature of H-bonding for 
the neutral and anion UQ species. Changes in H-bonding upon 
radical formation may be a mechanism for (de)stabilizing cofac- 
tors to fine-tune electron transfer processes in biological systems. 
Clearly it will be useful to undertake calculations similar to that 
presented here (modeling isotope edited FTIR difference spectra) 
for the UQ anion in the Qa binding site, and such calculations are 
underway in our lab. Based on the above one could argue that the 
methods used here, with only the quinone treated at the QM level, 
will be inadequate to simulate the experimental spectra associated 
with quinone anions in the Qa binding site. Or, for simulating the 
vibrational properties of quinone anions in the Qa binding site it 
will be necessary to treat key amino acids at a quantum mechan- 
ical level. Again, calculations are underway in our lab to test this 
proposal. 

For the neutral state of quinones occupying the Qa binding 
site, calculations with the H-bonding amino acids treated at the 
MM level lead to calculated spectra that are in excellent agree- 
ment with experimental spectra (Figures 3-5). Clearly, treating 
H-bonding amino acids quantum mechanically will not lead to 
improved modeling of the experimental spectra. We have in fact 
undertaken calculations in which neutral UQ and the H-bonding 
amino acids are treated using QM, and we have found that the cal- 
culated spectra are very similar to that obtained when only neutral 
UQ is treated using QM (and the surrounding amino acids are 
treated using MM) (Zhao et al., J. Phys. Chem. B in press). Thus, 
it is very clear that, at least for the case of neutral quinones occu- 
pying the Qa binding site, and as far as modeling isotope edited 
FTIR spectra is concerned, QM/MM calculations with only the 
quinone treated at the QM level need to be considered. 

An experimental ls O isotope edited FTIR difference spectrum 
for UQ in the Qa binding site is shown in Figure 5B. Three 
positive bands at 1660, 1629, and 1601 cm -1 are observed. By 
considering FTIR difference spectra obtained using PBRCs with 
unlabeled and specifically 13 Q and 13 C4 labeled UQ occupy- 
ing the Qa binding site, it was concluded that the 1660 and 
1601 cm -1 bands are due to the Ci=0 and C4=0 vibrations of 
unlabeled neutral UQ, respectively. It was also concluded that 
the 1628 cm -1 band is due to a UQ C=C vibration (Breton 
et al., 1994c). Since the C4=0 mode was so massively down- 
shifted (from ~1660cm _1 for UQ in solvent to 1601 cm -1 for 
UQ in the Qa binding site) it was suggested that this group must 
be engaged in very strong hydrogen bonding, presumably with 
HisM219 (Figure 2) (Breton et al., 1994c). This conclusion is 
difficult to rationalize based on the crystal structural data and 
other experimental data [see Wraight and Gunner (2009) for a 
review] . Such a conclusion is also not supported by the data pre- 
sented here. Specifically, the Ci=0 and C4=0 modes of DMNQ 
and VK are found at 1653 and 1643 cm -1 (Table 2), respectively, 
compared to ~1662cm~' for the coupled C=0 vibration in 
solution (Bandaranayake et al, 2006). Thus, the Ci=0/C4= O 
mode of DMNQ or VK in the Qa site is downshifted 9/19 cm -1 , 



respectively, compared to that found in solution. Such shifts sug- 
gest that both C=0 modes of DMNQ or VK are H-bonded in the 
Qa site, albeit quite weakly. 

From the experimental spectrum of VK (or DMNQ) in 
Figure 3B two positive bands are observed at 1651 and 1640 cm -1 
and one negative band at 1620 cm -1 . The two C=0 modes 
of unlabeled VK give rise to the positive bands at 1651 and 
1 640 cm~ 1 , but only a single band is observed at 1 620 cm~ 1 upon 
18 O labeling. Two interpretations for these observations have 
been proposed (Breton et al., 1994b). One suggestion is that upon 
18 O labeling the 1640 cm -1 band downshifts to ~1620cm _1 , 
while the 1651 cm -1 band downshifts to near 1640 cm -1 and 
decreases in intensity. The negative band near 1640 cm -1 due to a 
C = ls O group of VK is then masked by the positive band (also at 
1640 cm -1 ) due to the unlabeled C=0 group. A second hypoth- 
esis is that the two C=0 modes of unlabeled VK (at 1651 and 
1640 cm -1 ) both downshift to ~1620cm _1 upon ls O labeling. 
The different ls O induced shifts of the two C=0 modes results 
from differential coupling to C=C modes. 

The calculated data presented here allow us to address which 
of these interpretations could be appropriate, or if either is 
appropriate. The ONIOM calculations show that the Q=0 and 
C4=0 modes of unlabeled VK occur at 1652 and 1642 cm -1 , 
and that neither of these modes are coupled to C=C modes 
(Table 2). Upon ls O labeling the Ci=0/C 4 =0 mode downshifts 
24/18 cm -1 , respectively. The modes of ls O labeled VK also dis- 
play considerable coupling with C=C modes. Upon ls O labeling 
the Ci=0/C4=0 group couples with C=C r i ng in-phase/out-of- 
phase vibrations, respectively. Coupling of the C4=0 group to the 
out of phase C=C vibration leads to a large intensity enhance- 
ment, while coupling of the Ci=0 group to the in-phase C=C 
vibration leads to a large intensity decrease (Table 2). So, the cal- 
culations indicate that two separate uncoupled C=0 modes in 
unlabeled VK give rise to predominantly a single mixed mode 
(that carries most of the intensity) in ls O labeled VK. These calcu- 
lated results indicate that neither of the two previously proposed 
interpretations of the experimental spectra is correct. Clearly, the 
calculations presented here allow a more detailed insight into 
the nature of the bands in the experimental isotope edited FTIR 
difference spectra. 

Experimentally, the vibrational modes of DMNQ are at a 
slightly higher frequency (~3 cm -1 ) than corresponding modes 
of VK. Presumably replacing the C6 methyl group with a phytyl 
unit causes this difference. Interestingly, this small frequency dif- 
ference in the modes of DMNQ and VK is also found in our 
ONIOM calculated spectra (compare spectra in Figures 5C,B). 
This result is not entirely specific to the ONIOM method, how- 
ever, as a small shift is also found in the gas phase calculations 
(Figure 3D). 

The Ci=0 and C 4 =0 modes of DMNQ and VK are calcu- 
lated to be separated by 10 cm -1 . This separation cannot be due 
to differences in the molecular group attached at C&. It must be 
due to differences in how the two C=0 groups interact with the 
protein. Similarly, the two C=0 modes of DQ are calculated to 
be separated by 13 cm -1 , and this separation is also likely due to 
differences in how the two C=0 groups of DQ interact with the 
protein. 
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For UQi (and MQo) the separation of the C=0 modes is 
~32cm~ 1 (1660-1628). Some of the differences in frequency of 
the C=0 modes of UQ are due to the different orientation of the 
methoxy groups. If we assume that protein interactions with the 
C=0 groups gives rise to a 13 cm -1 separation in the frequencies 
of the two C=0 modes then this would indicate that the differ- 
ence in the orientation of the methoxy groups of UQ (or MQo) 
gives rise to a frequency difference of 19 cm -1 for the two C=0 
groups. This result is in quite good agreement with results from 
experimental spectra of UQ in solution, which show that the two 
C=0 modes are separated by ~16cm _1 (Breton et al., 1994c; 
Brudler et al., 1994). 

Normal mode vibrational frequencies are governed by molec- 
ular bonding force constants. These force constants relate to the 
electronic structure of the molecule. Since the calculated and 
experimental spectra for DMNQ and VK are virtually the same 
it is concluded that the replacing the phytyl unit at C6 of VK with 
a methyl group does not appreciably perturb the electronic struc- 
ture of the NQ ring. In addition, the tail does not perturb the 
protein in a way that significantly modifies any pigment protein 
interactions. 

In VK and UQi, the "kink" in the hydrocarbon chain after 
the first carbon atom is 3.1-3.2° higher in ONIOM calcula- 
tions compared to gas phase calculations. It seems therefore, that 
the hydrocarbon chain is somewhat constrained relative to the 
quinone ring when incorporated into the Qa binding site. It is 
not clear if this is a significant constraint, however. The Ci=0 and 
C 4 =0 bonds of MQ 0 and UQi (and of DMNQ and VK) are vir- 
tually unaltered in ONIOM calculations (Table 1). The hydrogen 
bond lengths for the C=0 groups are also little altered (Table 1). 
The distance of the C4=0 oxygen to the non heme iron atom 
is 0.029/0.055 A longer for UQi/VK compared to MQn/DMNQ, 
respectively, suggesting a very small change in orientation of the 
UQi/VK head-group (since the iron atom is fixed) compared to 
MQn/DMNQ in the Qa binding site. The hydrocarbon chain may 
therefore, lead to a very slight change in the orientation of the 
quinone ring in the Qa binding site. Figure 2F shows that there is 
only a very small difference in the orientation of DMNQ relative 
to UQi. 

The experimental ls O isotope edited spectra for MQo and 
UQi are quite different (compare Figures 5B,C). This difference 
cannot be modeled computationally (Figure 5A). The experi- 
mental isotope edited spectrum for MQo is considerably noisier 
than the spectrum for UQ (Breton et al., 1994a). As far as we are 
aware the Qa /Qa FTIR difference spectrum for RCs with MQo 
in the binding site have never been reproduced, so the accuracy of 
the spectrum may be somewhat questionable. On the other hand 
MQo may be able to incorporate into the Qa binding site with the 
methoxy groups oriented in several different ways. As described 
previously, each of these methoxy group conformers will have 
slightly different spectra (Lamichhane et al., 2010). This hetero- 
geneity in orientation of MQo in the Qa binding site may be a 
factor that contributes to differences in the experimental spec- 
tra for MQo and UQ, as demonstrated in Figures 5B,C. In spite 
of this, there are some overall similarities in the MQo and UQi 
experimental spectra in Figures 5B,C. For the MQo spectrum 
positive bands are found at 1665, 1631, and 1608 cm -1 . For UQi 



the bands are at 1660, 1628, and 1601 em -1 . For the MQ 0 spec- 
trum negative bands are found at 1616 and 1602 cm -1 . For UQ 
the bands are at 1613 and 1586 cm -1 . 

The ONIOM calculated and experimental isotope edited spec- 
tra for VK are very similar, although the intensity ratios of the 
different bands do not appear to match well. For example, the 
intensity ratio of 1652 and 1644 cm -1 bands in the calculated 
spectrum (Figure 3 A, dotted) is ~2.1, compared to a ratio of 
~0.95 in the experimental spectrum. To investigate this further 
we have also calculated the 13 C isotope edited spectrum for VK 
in the Qa binding site (Figure 6A) (global 13 C labeling of VK), 
and compared it to the corresponding experimental spectrum 
(Figure 6B). For completeness the calculated 13 C isotope edited 
spectrum for VK in the gas phase is shown in Figure 6C. Again, 
it is evident that the calculated gas phase spectrum in no way 
resembles the experimental spectrum (compare Figures 6B,C). 
The normal modes (frequencies, intensities, and PEDs) that give 
rise to the various bands in the ONIOM calculated 13 C isotope 
edited spectra of VK in Figure 6A are listed in Table 3. 

Notice that the ONIOM calculation predicts that the two 
C=0 modes remain completely separate (do not mix) upon 
13 C labeling (Table 3), which is unlike the behavior observed 
for the C=0 modes upon ls O labeling (see above). Also notice 
that in the experimental isotope edited spectrum (Figure 6B) 
the 1640 cm -1 band is now considerably more intense than 
the 1651 cm -1 band, whereas in the calculated spectrum the 
higher frequency band is considerably more intense. The ori- 
gin of these mode intensity differences are not entirely clear and 
are currently being investigated by considering calculations asso- 
ciated with VK in the presence of various types of H-bonding 
molecules. 

CONCLUSIONS 

We have shown that ONIOM type QM/MM calculations can 
be used to model experimental isotope edited FTIR difference 
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FIGURE 6 | (A) ONIOM calculated 13 C isotope edited DDS for neutral VK 
in the Q A binding site. Experimental spectrum is also shown (B), taken 
from Breton (1997) (with permission). (C) DFT calculated 13 C isotope edited 
DDS for neutral VK in the gas phase. ONIOM and gas phase calculated 
spectra were scaled by 0.9718 and 0.9608, respectively. 
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Table 3 | Normal mode frequencies (in cm 1 ), intensities (in km/mol) and PEDs (in %) calculated using ONIOM methods for unlabeled and 13 C 
labeled neutral VK. 
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Frequency shifts upon ,3 C labeling are also listed. Negative signs in the PEDs refer to the relative phase of vibration of the internal coordinates. Only internal 
coordinates that contribute at least 5% are shown. Mode frequencies were scaled by 0.9718. 



spectra obtained using PBRCs that have had several different 
quinones incorporated into the Qa binding site. The fact that the 
many different spectra can all be modeled is a clear indicator of 
the appropriateness of the approach. 

The calculated spectra appear not to depend on whether the 
quinone incorporated has a prenyl/phytyl unit or a methyl group 
attached at C$. The electronic structure of the quinone ring is 
therefore, not sensitive to the presence or absence of a hydro- 
carbon side chain at This suggests that the hydrocarbon side 
chain does not significantly constrain the quinone ring in the Qa 
binding site. 

Comparison of the calculated and experimental spectra, in 
combination with a consideration of the calculated PEDs of the 



normal modes, allows a direct assessment of the appropriate- 
ness of previous suggestions as to the origin of the bands in the 
experimental spectra. 
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